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ABSTRACT 

The flux-injection hypothesis for driving coronal mass ejections (CMEs) requires the trans- 
port of substantial magnetic energy and helicity flux through the photosphere concomitant with 
the eruption. Under the magnetohydrodynamics approximation, these fluxes are produced by 
twisting magnetic field and/or fiux emergence in the photosphere. A CME trajectory, observed 
2000 September 12 and fitted with a flux-rope model constrains energy and helicity budgets for 
testing the flux-injection hypothesis. Optimal velocity profiles for several driving scenarios are es- 
timated by minimizing the photospheric plasma velocities for a cylindrically symmetric flux-rope 
magnetic field subject to the flux budgets required by the flux-rope model. Ideal flux injection, 
involving only flux emergence, requires hypersonic upflows in excess of the solar escape velocity 
617kms~^ over an area of 6 x 10® km^ to satisfy the energy and helicity budgets of the flux-rope 
model. These estimates are compared with magnetic field and Doppler measurements from Solar 
Heliospheric Observatory /Michelson Doppler Imager on 2000 September 12 at the footpoints of 
the CME. The observed Doppler signatures are insufficient to account for the required energy 
and helicity budgets of the flux-injection hypothesis. 

Subject headings: Sun: coronal mass ejections - Sun: photosphere - Sun: surface magnetism 



INTRODUCTION 



Over the last 10 years, the flux-rope model developed by [Chen (1989 19961 has been used to describe 
the dynamics of coronal mass ejections (CMEs) observed by the Large Angle Spectrometric Coronagraphs 
(LASCO) aboard the Solar and Heliospheric Observator-^{SOHO) (Chen et al. 



1997 



2000 



Wood et al. 



1999 Krall et al. 2001 Chen et al. 2006 Krall et al. 2006[ ). However, the photospheric fiux injection paradigm 
used to initiate and drive the eruption has been criticized because the surge of electromagnetic energy fiowing 
through the photosphere is "difficult to reconcile with the extremely tranquil conditions that exist during 



fiares and CMEs" (Forbes 2000 2001 Rust 2001). Chen et al. (2000), Chen (2001), Krall et al. (20011, 



Chen & Krall (2003) and Chen & Kunkel (2010) have attempted to address these criticisms and a previous 



study has examined the implications of uniformly twisting to coronal footpoints of the flux rope ( Krall et al 
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2000). However, no quantitative comparisons between flux- injection hypothesis and detailed photospheric 



observations have been considered by the formal literature. 

The goal of this paper is to provide a framework for testing the flux-injection hypothesis through the 
photospheric signatures implied by the energy and helicity budgets of CMEs described by the flux-rope 
model. The paper is organized as follows: Section [2] describes the flux-rope model of Chen ( 1989 ) and a 
simple extension of the flux-rope model magnetic field into the photosphere. 

Section |3] develops the photospheric fluxes necessary to satisfy the energy and helicity budgets of CMEs 
fitted with the flux-rope model. The photospheric magnetic field combined with the photospheric fluxes is 
used to estimate minimum velocities necessary to satisfy the energy and helicity budgets required by CME 
trajectories fitted by the flux-rope model under the flux-injection hypothesis. Two examples of CME trajec- 
tories fitted with the flux- rope model are used to constrain the photospheric velocities. 1) The first event is 
the 2000 September 12 CME that erupted from decaying NOAA active region 9163. The height-time data 



for the CME trajectory (Chen et al. 2006) were derived from measurements of the filament in absorption 



observed by the Global Ha Network at Kanzehohe Solar Observatory (KSO) in Austria (Steinegger et al 



2000), and LASCO C2 and C3 observations of the filament in emission ( Brueckner et al.|1995 ). Complemen 



tary observations of the filament were made in Fe XII 195 A by the SOHO /WY Imaging Telescope (EIT) 
instrument ( Delaboudiniere et al.|199"5 ). The eruption occurred shortly after 11:30 UT and the filament first 
appeared in LASCO C2 at 12:30 UT. This CME was associated with an Ml.O class fiare with distinct flare 
ribbons that persisted for 2 hr. Various aspects of this event, such as morphology, timing, and reconnection 



rate have been discussed by Vrsnak et al. (20031, Schuck et al. (2004b) and Qiu et al. (2004). 2) The second 



event is the 2003 October 28 CME which originated from large complex NOAA active region 10468 at about 
11:00 UT. The event was extremely fast, and thus the height-time data for the leading edge of the CME 
consists of only one LASCO C2 image at about 11:30 UT and four subsequent LASCO C3 images. This event 
was associated with an extremely powerful X17 fiare and consequently has received extensive attention in 



the literature ( 


Seppala et al.|2004 


Skoug et al.|2004 


Woods et al.||2004 


Zurbuchen et al.|2004 


Bieber et al. 


|2005||Chi et al. 2005 ; 'Degenstein et al.||2005| |Hu et al. 2005 


Gopalswamy et al. 2005; Looper et al.||2005| 



|Panamraj u~^Chakrabarti|2005[[Tsurutani et al.|2005[|Yurchyshyn et al.|2005[[KraU et al.|2006[ [Manchester] 
[et al.,2008] ). 

Section [4] compares the photospheric velocities implied by the flux-rope CME trajectory event from 
2000 September 12, against detailed photospheric Doppler measurements from the Michelson Doppler Imager 
(MDI) aboard SOHO (Scherrer et al. 1995 1. Finally, Section [s] compares these theoretical results and Doppler 
observations with previous work. 



THE FLUX-ROPE MODEL 



Figure [T] shows a schematic diagram of the flux-rope CME model current loop, adapted from |Chen| 
(1989 1996[ ). [Shafranov (1966) derived the magnetohydrodynamic (MHD) forces per-unit-length acting the 
major TZ and minor a radii of a current carrying toroidal section 



c2 7^ 



In 



?7^ 

a 



/a — 



2 

({B^ 
Bl 



2Bl 



2 7^ Be 

a Boa 



- 1 



-l + A/3 



(la) 
(lb) 



- 3 - 



Be 




Fig. 1. — Schematic diagram of the flux-rope model current loop adapted from Chen (1989 1996). The 



subscripts "(/)" and "0" refer to the toroidal and poloidal directions in the cylindrical {R, 4>, Z) and the local 
polar (r, 6, Q coordinate systems for a circular torus where Z is out of the plane of the page. 



where is the toroidal current, c is the speed of light, (-B^) is the average toroidal field inside the current 
channel, is a prescribed function describing the overlying coronal field perpendicular to the flux rope, 
Bea = Bg (a) = 2I^/ac is the poloidal magnetic field at the edge of the current channel, and 



2j;drrBUr) 



c 

2P 



4> Jo 



drrBl {r)^0{l), 



(2) 



is the internal inductance. The quantity A/3 = 8tt{{p) — Pc)/Bg^ is the differential of plasma f3 based on 
Bga where (p) is the average internal flux-rope pressure and Pc is the coronal pressure. The subscripts "(/)" 
and "0" refer to the toroidal and poloidal directions in the cylindrical {R,(j),Z) and the local polar {r,9X) 
coordinate systems for a circular torus with R = TZ + r cos 6, Z = r sinO and C = <t> where Z is out of the 
plane of the page. Chen ( 1989 ) applied Equation ([!]) at the peak of the flux rope to derive the dynamical 
equations 
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where iJ > is the height of the center of the current channel above the photosphere henceforth referred 
to as the "apex." The *'s in Equation (|3| denote quantities evaluated over the cross section at the apex 
where M = (n*) mi: is the mass per- unit-length of the flux rope and (n*) is the average density of the 
flux rope. The additional terms 



Id ^CdUcm a* (Vsw - V) I Vsw - V\ 



(4a) 



and 



/g = Tralmg {uc - (n*)) , 



(4b) 



introduced in Chen ( 1989| ), represent the drag and gravitational forces per- unit-length of the flux rope 
respectively where Cd ~ C (1) is the drag coefficient, rn ~ 1.67 x 10^^"* g is the mass of hydrogen, ric is the 
ambient coronal density, V = dH/dt, Vsw is the velocity of the ambient solar wind, and the gravitational 
acceleration is 

ff- 2' (5) 

with = 2.74 X 10'* cms~^ and Rq = 6.74 x 10^" cm. The flux-rope footpoints are separated by a distance 
Si in the photosphere and assumed to remain fixed throughout the evolution of the flux-rope CME by the 



"dense subphotospheric plasma" (Chen 1989). The current-channel radius at the base of the corona Uc is 



also assumed to remain fixed throughout the evolution of the flux rope. The major radius TZ and height H 
above the photosphere are related by 



n EE ^' and H >0. 

2 H 

The length of the flux rope above the photosphere is L = 2 tt 72. where 

Q ^ / 1 - ^/^ H> St/2, 



ip/lT 



H < Si/2, 



(6) 



(7) 



and if = arcsin {Si/2H) ( |Chen||1989] |Krall et al.||2000[ ). Note the variables that describe the plasma at the 
apex of the flux rope (n* ) , (p* ) , Bga, i i and will evolve with time and parameters of the interplanetary 



medium ric. Be, and Vsw are implicit functions of H. (see Chen 1996 for a complete description of the 



interelations between variables and a model for the parameters of the interplanetary medium) (Krall et al. 



2000) 



Flux-rope equilibria are determined from 



0. 



(8) 



Chen & Garren (1994), Chen (19961, Krall et al. (2000), and KraU & Chen (20051 proposed local polar 



equilibrium profiles for 5^ and large aspect ratio TZ/a » 1 

Be (r) = Bea 
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compatible with the flux-rope model with a toroidal flux of 
and carrying a bare toroidal current of 



(10a) 



I^ = I^ = B^aac/2. (10b) 
There is no return current which is consistent with the large-scale current systems in active regions analyzed 



by Wheatland (2000). Since Equations (9a I and (9b I do not represent a force- free equilibria JxB = 0, the 



poloidal and toroidal magnetic flelds and corresponding currents are decoupledj^ Consequently, the cutoff 
of the toroidal flcld i?,^ at r = a is arbitrary relative to the poloidal magnetic field. For example, 



B^ (r) = B^a e 



2 / 2 

— r /a 



(11) 



is a valid toroidal field model and contains the identical amount of flux as Equation ( 9b ) 



The toroidal energy per-unit-length of the flux rope is 

Uc = - ( drrBfir)- —Bf^a'^, 



(12a) 



and the poloidal energy per unit length contained within a distance of radius r of the toroidal axis of the 
flux rope is 

Ue (r) = \f^ dr' r' B^ (r') ^ [73 + 120 log (rja)] . (12b) 



The local polar equilibrium profile Equations ( 9a I and ( 9b ) does not formally admit a bounded poloidal 



energy per unit length along the flux rope — a necessary physical condition for admissibility — because Bg ~ 
and Ug (r) ^ log (r) for r > a. However, for the closed circuit representative of the schematic flux rope 
sh own in Figure [T] the inte gral is cutoff at distances of order of the dimension of the circuit (see pp. 136-141 
in Landau fc Lifshitz|1960 ) which is roughly half the footpoint separation r^. ~ S'f/2. The fraction of poloidal 



magnetic energy contained inside radius r is then 



6e ir,a,rc) 



_ Ug{r) 73 + 120 log (r/a) 



Ug (re) 73 + 120 log [rja) ' 



(13) 



Chen et al. (2000) have argued that a "large fraction of the injected poloidal energy is in the magnetic field 
outside the current channel," denoted by r = a. Indeed, in the corona, roughly 1/3 of poloidal energy is 
contained within the current channel r < a. However, more than 2/3 of the poloidal energy is contained 
within the region r < 2 a near the current channel. In the photosphere, the amounts are reduced because the 
current channel is narrower, but the poloidal energy is contained within the region r < 2 a near where the 
current channel remains significant varying from 1/3 2/3 depending on the strength of the toroidal field 
in the photosphere. The comparison of the flux-injection hypothesis against photospheric observations in 
Section [4] is limited to the region r <2a where the poloidal energy transport is substantial and the poloidal 
magnetic field in significant. 



Although the magnetic field model Equations ( 9a ) and ( 9b ) was originally proposed as an approximate 



local description for large-aspect ratio coronal toroidal flux-rope magnetic fields, it has also been used as 
an initial equilibrium for investigating photospheric signatures in simulations of subphotospheric fiux ropes 



^In equilibrium, the magnetic pressure is in detailed balance with gravity and kinetic pressure. 
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Fig. 2. — Schematic diagram of the flux-rope leg through the photosphere adapted from Chen (2001 ). The 
current channel carrying J^, the top of which is shaded, and the poloidal field Bg, and two community field 
lines extending into the corona Be are shown. The subscripts "c" and "p" refer to coronal and photospheric 
values respectively. 



( Chen fc Huba|2005a|b 2006 ) and similar models have been implemented for modeling the photosphere flux- 
rope dynamics ( Sakai et al.|2000 Sakai et al.|2001 ). Figure [2] shows a schematic diagram of the flux-rope leg 
through the photosphere adapted from Chen (2001 ) in the local polar r, 6, ( coordinate system. The vertical 



current channel carrying Iq, the top of which is shaded, and the private poloidal field Bg, and two community 
field lines extending into the corona Be are shown. The subscripts "c" and "p" refer to values at the base of 
the corona and in the photosphere, respectively. Two general constraints may be applied quasi-adiabatically 
to extend the coronal model into the photosphere: conservation of vertical flux and vertical current ( |Krall] 
fc Chen|[2005l ) 



B^p 
Bgp 



B. 



Bo 
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(14a) 
(14b) 



Combining these relationships with Equations ( 9a ) and ( 9b I produces an approximate local model for the 
magnetic footpoints of the flux rope in the photosphere with _B^c = B^c at the base of the corona. The details 
of the flux-tube expansion in the chromosphere is not explicitly specified in this approximation. However, 
the expansion implies a local radial field Br (r) to balance the decrease of B^ {() with height under the local 
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constraint of V • B = 0. 



THE PHOTOSPHERIC FLUX BUDGETS OF THE FLUX-ROPE MODEL 



The flux-rope model requires specification of two quantities in addition to the geometrical aspects of 
the flux rope, namely the toroidal and poloidal fluxes ( Chen||l996 ) 



$0 = TJ" flc = '^oiistant, 



where 



£=i^<log(8R) 



^4 



[a* log (a*) - flc log(ac)] 



(15a) 
(15b) 

(16) 



is the inductance (Landau & Lifshitz 1960[ Krall et al. 2000). The toroidal energy is estimated from the 
arclength of the flux rope above the photosphere}^ 



C/0 



97re 7^ 



r2 „3 



9 n 



constant, 



20 a* " 20n Qc a 
whereas the poloidal energy is related to the inductance and poloidal flux 



The toroidal flux 



is conserved because the toroidal field i?^ 



(17a) 
(17b) 

and current-channel radius 



flc at the base of the corona are held constant in time during the eruption — B^a and the minor radius a 
may vary along the flux rope subject to flux conservation. The toroidal energy Uij, is conserved because the 
flux rope erupts self-similar ly TZ/a^, w constant. Consequently, the quantities that will manifest dynamical 
changes in the photosphere during the eruption are the poloidal power dUg/dt and the rate-of-change in 
poloidal flux d$e / dt which are related to the photospheric energy and helicity fluxes respectively. 



The MHD induction equation derived by combining Faraday's Law 

dtB = -cVxE, 

with Ohm's Law 



to obtain 



a c 



dtB = Vx (vxB) -cVx 



(18) 
(19) 
(20) 



where v is the plasma velocity and cr is a spatially variable conductivity. The magnetic energy in the corona 



is formulated by dotting the induction Equation (20) with the magnetic field B and using Ampere's Law 
without displacement currents 



VxB 



J, 



(21) 



^The flux-rope minor current-channel radius varies linearly with arc length a{£) = Uc + (a» — ac) i from the footpoint at 
£ = to the apex at £ = 1. 
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to derive Poynting theorem 

dUu _ d r dV 
dt ~ dt Jv^^ 



S2 



1 

4 7r 



dSn 



Bx (vxB) + - (JxB) 



dV 



■ {JxB) 

c 



(22) 



where the volume integrals Vc are over the corona, chromosphere, and transition region, the surface integrals 



Sp + 5*0 



are over the photosphere Sp and the outer corona at R ^ R, 



'0y 



bounding the volume S 

and n is the surface normal pointing into the coronal volume (radially outward at the photosphere Sp and 
radially inward at the outer corona ^oc)- Assessing the energy budget of the region between the photosphere 
and outer corona by tracking the Poynting flux through the photosphere and the energy leaving the corona 
through eruptive phenomena provides an estimate for the free energy available for producing flares and 
CMEs (see Figure 1 in Kusano et al.|2002 ). The flrst and second terms in Equation (22 1 represent the ExB 
Poynting flux through the surfaces and the third and fourth terms represent the conversion of magnetic 
energy to kinetic energy though work done by the JxB force on the plasma and Ohmic heating through 
resistivity respectively. 



The hehcity for the flux rope may be writteij^ ( |Chen||l996| |Chen fc KraII||2003[ ) 



K : 



(23) 



Comparisons between integrated photospheric hclicity flux and the poloidal flux injection proflle for the 
flux-rope model can be made. However, because the field lines of the flux rope penetrate the photosphere, a 



gauge- invariant relative hclicity must be used for estimating hclicity fluxes through the photosphere ( Berger 
fc Field|[T984l ) 

AK^ I dV {A-B-Ar- Br) , (24) 

where Vc corresponds to the volume above the photosphere and J5r = V X Ar are the reference fields which 
are chosen to match the normal components of B and the tangential components of A respectively at the 
surface: 

{A-AR)xn\s = 0, (25a) 
(B - Br) • n\g = 0. (25b) 

These boundary conditions are sufficient for the equivalence of the relative helicity defined by [Berger ^ 



Field 


1984 


(1985 





AK 



[ dV (A 

Jv,. 



Ar) ■ {B ~ Br) , 



(26) 



because 



/ dV {Ar-B-A-Br) 



[ dV [Ar-VxA- A-VxAr) . 

I dVV-{AxAR), 

- / dSn ■ {AxAr) , 
Js 



0, 



(27) 



"See pp. 21,516 in 



Chen 



1 1996 I. This expression is exact in the large aspect ratio limit Ti/a 2> 1. 
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with Equation (25 I 



A judicious choice for the reference field is a potential field 
in the Coulomb gauge 



V • Ar = 0, 



with the additional boundary condition 



n ■ A 



For these conditions, this reference field has zero helicity 



Kn = / dVAn-Bn = / dVAn-V-^R = 



dSn-AR^R-l dV^RV-AR=0. 

Vc 



(28a) 

(28b) 
(28c) 

(29) 



Using Equation (24) or (26) with Equations (28a) and (28b) and boundary conditions (25) and (28c 



the Poynting theorem for the magnetic helicity into the corona then takes a particularly simple form 
dAK d 



dt 



dt 



dVA-B, 



2 (j) dSn- 



ArX {vxB) ArXJ 

a 



2c I dV 



J B 



(30a) 
(30b) 



The first and second terms in Equation ( 30b ) represent the helicity flux through the photosphere and the 



third term represents helicity dissipation in the coronal volume. 

Although the photosphere is not ideal, with a magnetic Reynolds number of Rm — U L/rj ~ fO^ — fO^ 
where U is the typical velocity, L is the typical gradient scale, and ry = c^/ {Ana) is the magnetic diffusivity, 
the ideal approximation has been demonstrated to be adequate for inferring plasma velocities from magnetic 
field dynamics in convection zone simulations]^ with magnetic Reynolds numbers as low as Rm ^ 10'^ ( Welsch 
et al.||2007[|Schuck|[2008| . The ideal MHD induction equation becomes 



dtB = Vx {vxB). 



(31) 



Within the flux injection paradigm, consistent with Equation (f5a), the toroidal magnetic field in the pho- 



tosphere does not change during the eruption ( Chen et al.||1997[ ) 

dtBi; = C • Vx (vxB) « 0. 



(32) 



This relationship also represents an observational constraint on the vertical magnetic field of the fiux rope in 
the photosphere because large changes in line-of-sight magnetograms near disk center have not been observed 
during eruptions. This constraint on the vertical magnetic field implies 

vxB^ Vhip + (vr Eg -~BrVe) C, (33a) 

where the subscript "h" refers to the horizontal (r, 9) coordinates of the local polar coordinate system and 



tj) is the electrostatic potential. Solving Equation (33a) for v produces 



V = s ^ h v\\ B \B\ 



(33b) 



= See 



Abbett et al. 



I 2004 1 
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where wy is the field- aligned plasma velocity. For a locally cylindrical flux rope with dg = and Br = 
consistent with Equation ([9| , the evolution of the radial and poloidal fields at the photosphere are determined 

by 



dtBr - 0, 

dtBg = drci' - B^^ (dcB^) (dr^j) , 

dtB^ = 0. 



(34a) 
(34b) 
(34c) 



Chen (1996) and Chen & Kunkel (2010) have noted that the poloidal magnetic field at the base of the 



corona may only increase by 20%-50% during the first 30 minutes of the eruptions and Krall et al. (2001) 
have shown an event that requires no perceptible increase in the poloidal magnetic field (see Figure 8(a) in 
Krall et al.||2001 1. To the extent that d(; « 0, these results are consistent with MHD. 



3.1. The Ideal Poynting Fluxes in the Photosphere 



From Equations (22 1, (30b), and (33b), the ideal MHD energy flux through the photosphere is 

1 



{vt; Bh - Bq Vh) ■ Bh, 



1 

Ait 



{CxVh^)-Bh, 



and the hclicity flux is 



k(; = 2 {vt^Bh- Bt^Vh) ■ Ar, 



(35) 



(36) 



Equations ( 33 ) , ( 34b ) , ( 35 ) , and ( 36 ) are general with respect to the constraint ( 32 ) ; the fluxes through the 
photosphere are completely specified by the electric potential tp the magnetic field. The flux-injection 
hypothesis has always appealed to some "unspecified sub-photospheric process" for initiating and driving 



the CME (Krall et al.| 


2000| see also Chen"1989t 'Chen & Garren||1993| 


Chen||1996! 'Chen et al.||1997 


Chen 


1997 


Chen et al.|2000 


Chen|2001 


Chen & Krall^2003j Krall et al.|2001 


Chen & Kunkel 20101. Nonetheless, 



in the ideal MHD limit, energy is transported through the photosphere by the term ^ • Bx (vxB) and 
helicity is transported by C • ArX (vxB). These terms integrated over the footpoints in the photosphere 
(including the region outside the current channel) must balance the poloidal power requirements dUg/dt 
and helicity requirements dK/dt ~ $0 d^g/dt of a CME described by the flux-rope model. Temporally, the 
increase in the fluxes at the footpoints should precede the eruptions by at least (5i ~ 4 minutes to account 
for the transport of magnetic field along the flux rope to the apex of the CME at the coronal Alfven speecj^ 
Va ^ 1-5 X lO^cms-i. 

Assuming azimuthal symmetry dg = 0, the poloidal power injected through the flux-rope footpoints can 



be written in terms of Equation (17b) or (35) 



dUg 
dt 



2 



dt 



1 dC\ $2 



-2 X 



dr r Bg drtp ~ 



dr r Bg drip, 



(37) 



®The Alfven speed is based on fiux-ropc initial conditions from (Chen 19961 with (n) ~ 7.5 X lO'^ cm ^ and _B ~ 6 G. Note 
that this estimate is an order of magnitude larger than estimates corresponding to the same height range by [Regnier et al.| 
l 2008| . 



and because the toroidal flux is constant, the rate of change of helicity is 



dK 

~dt 



dt 



-2 X An j dr r A^ie dr'4' y 
'o 



(38) 



where the fluxes are assumed to be equipartitioned between the two footpoints. Equations (37) and (38) 



represent independent constraints between the flux-rope dynamics and the energy and helicity fluxes in the 
photosphere. 

Balancing the photospheric Poynting fluxes against the poloidal energy and helicity budgets of CME 
trajectories fit by the flux-rope model requires choosing a photospheric model for the magnetic field and a 
velocity profile or drip- For the former, the simple photospheric extension of the flux-rope model described 
in Section [2] is used. For the latter, the optimal velocities for transporting energy and helicity across 
the photosphere consistent with the respective overall budgets of the flux-rope model are estimated via 
constrained variational calculus (see Appendices |A] and [b]). Optimal, in this context, means the minimum 
root-mean-squared (rms) photospheric velocities corresponding to the minimum photospheric kinetic energy 
(for constant density). 

There are two physical mechanisms for transporting energy and helicity through the photosphere and 
into the corona: 



1. Twisting the magnetic field in the photosphere through poloidal motion. Under this mechanism and 
the magnetic field model in Equation all the energy and helicity must be transported through 
the photosphere and into the corona within the current channel where ^ 0. Minimizing J drrvg 
with = 0, the rms poloidal photospheric velocity inside r < ap constrained by the energy budget is 
derived in Appendix \K\ 



„2\l/2_ 

and constrained by the helicity budget derived in Appendix [B] 
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4V70 
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CP I 
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2\l/2 



d^e 
dt 



1 



70 1 

437 v/|i?CcScpl 

The twisting of footpoints does not change the net current carried by the flux rope (see discussion on 



pp. 



30 and Appendix C| 



Flux injection involving the emergence of poloidal flux transported through the photosphere by vertical 
motion. Under this mechanism and the magnetic field model in Equation ([£]), the energy and helicity 
transport are not limited to the current channel, but significant transport occurs within about r < 2 ap. 
Minimizing J drr with vg = 0, the rms vertical photospheric velocity inside r < 2ap constrained by 
the energy budget is derived in Appendix [A] 



2a„ 



dUg 



dt 



V12516735 \B^p\ 



2Bi 



2998 



Cp' 



1155 B, 



Cc' 



A13 



Note that Equation (A13) approaches a constant asymptotically if either Tc — ?> cxd or — )• oo, while 
holding the coronal flux constant i?^c 
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and this value is entirely determined by the current-channel radius and the poloidal magnetic field at 
the base of the corona. For efficient energy transport, a significant fraction of the energy transport for 
flux emergence will occur in and near the current channel. The corresponding rms vertical photospheric 
velocity inside r < 2 Op constrained by the helicity budget is derived in Appendix [B] 



2\l/2 



/ 2\-^/^ 
W2a 



dt 



V12516735 \B^p\ 



8 Bsc Oc 1 2998 B(;pr^ - 1155 B^;^ a? 



For comparison, minimizing J drr (drVQ)"^ with vg = produces the minimum shear estimate (constant 
velocity r < fc) which has rms vertical photospheric (everywhere) of 



/ 2\l/2 



dUe 



for the energy constraint and 

.2\l/2 



dt 



d<^fi 



120 



al 73 + 120 log (B^p rjB^c Qc) ' 



dt 



30 



Bee ac 73 + 120 log (B^p re/B^c flc) ' 



(40a) 



(40b) 



for the helicity constraint. 



3. Both limits above require wasted flows along field lines f || ^ 0. However by minimizing J drr [vg + w^) 
with wy = 0, the most efficient rms total velocity inside r < 2 Cp constrained by the energy budget is 
derived in Appendix [A| 



\ I la 



dUe 



dt 



2 Bt^p rl V23IO v/21674B2^ + 14421 5^;^ B, 



Cp 



i2 \Bec B^p (23984 Bl + 14421 B^^ ^Cp) ^? " 9240 a* B^ B, 



A17a 



and the most efficient rms total velocity inside r < 2 ap constrained by the helicity budget is derived 
in Appendix [B] 



/ 2\l/2 



dt 



B^p rl ^/2310 v/21674 B^^ + 14421 B^^ B^p 
2 oc \B^p (23984 Bl + 14421 B^^ B^p) rl - 9240 Bl^ B^d ' 



Bll 



where vn = everywhere. 



These three driving scenarios are used to estimate the expected rms velocities near the current channel of 
CME flux-rope footpoints. 



3.2. Velocity Estimates Constrained by CME Trajectories 



The flux-rope model has been fitted to the CME trajectory observed on 2000 September 12 (Chen 



et al. 2006 Chen & Kunkel 2010). 



Chen et al. (|2006|) and Chen & Kunkel (|2010|) find that the observed 

^ and an energy 



dynamics is consistent with a peak poloidal flux injection rate of d^^g/dt = 1.4 x 10^^ Mxs 
injection of 2 — 4 x 10"^^ erg over 40 minutes (James Chen, personal communication 2008 November). The 
average power requirement of the flux-rope model (dUg/dt) ~ 1.3 x lO^^ergs"^ is the lower limit for the 
power budget of the CME over 40 minutes — the instantaneous power budget is likely to exceed this value 
dUg/dt > (dUg/dt). The separation distance between the filament footpoints and the flux-rope footpoints 
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Fig. 3. — left: photosphcric current-channel radius Op (solid line) and poloidal magnetic field Bgp (dashed 
line) as a function of the average vertical magnetic field in the current channel B,^p assuming coronal param- 
eters consistent with the 2000 September 12 CME trajectory: Oc = 7.5 x 10® cm, Bee — Bqc — 4 G. Right: 
photospheric magnetic field profiles for B^p = 500 G 



were estimated to be Sf : 
channel radius Oc = {S{ 



3.5 X 10^° cm and S'f ~ 5 x lO^'^ cm respectively which leads to a coronal current- 
- Sf) /2 ~ 7.5 X 10® cm. The magnetic field at the base of the corona is estimated 



to be Bsc = Bqc — 4 G ( Chen fc Kunkel|[2010 ). Figure |3] shows the photospheric current-channel radius Op 
(solid line) and poloidal magnetic field Bgp (dashed line) as a function of the average vertical magnetic field 
in the current channel determined from Equations (14a I and (14b). 



Figure |4] shows velocity estimates for the 2000 September 12 CME trajectory modeled by |Chen et al. 



(2006) and [Chen fc Kunkel| (2010) as a function of the average vertical magnetic field in the current channel 
using a cutoff scale Tc = S'f /2 = 2.5 x 10"'^° cm. The magnitude of the vertical magnetic field can be interpreted 
as height with B^p = = 4 G corresponding to the base of the corona and B^p ~ 500 G corresponding 
to the photosphere. The solid lines correspond to rms velocities estimated from the energy budget and the 
dashed lines correspond to rms velocities estimated from the helicity budget. In the left panel, the blue and 
red lines are the optimal twisting (vgYj'^ and emergence velocities at the flux-rope footpoint inside 

r < Op and r < 2 Cp respectively, corresponding to the minimum velocity solutions in Appendices [A] and |B] 
The green lines indicate the emergence velocity corresponding to the minimum shear (constant velocity 

1/2 

r < Tc) solution. In the right panel, the blue, red and black lines correspond to the poloidal {vf)^^ and 

vertical {vl)]^"^ components of the optimal total velocity respectively with U|| = at the flux-rope 

footpoint inside r < 2ap. In both panels, the dotted black line is the solar escape velocity Vq — 617kms~"'^, 
the dotted green line is the photospheric sound speed Cg = \J 1V\>I P\i — 7.2 km s^^ with 7 
density pp ~ 5.85 x 10~* g/cm'^, and pressure pp 
sun (Vernazza et al. 1981) interpolated to the r 



MDI (Jones 



1989 



inside r < za 



5/3, mass 

1.82 X 10^^ gcm~^ s^^ from VAL-C model for the quiet 
_______ = 1 height of 240 km for Ni I 6767.8 A hue imaged by 

Bruls|l993 ). The dotted blue hue is the average Alfven speed (VA)2ap = (^)2ap / ^74 tt pp 



he velocities based on the energy and helicity budgets are in close agreement and this 
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Fig. 4. — Optimal velocity estimates for the 2000 September 12 CME trajectory modeled by |Chen et~al 



(2006) and Chen & Kunkel (2010) as a function of the average vertical magnetic field in the current channel. 



The solid lines correspond to rms velocities estimated from the energy budget and the dashed lines correspond 
to rms velocities estimated from the helicity budget. Left: the blue and red lines are the optimal twisting 
(vg)]/'^ and emergence {v^f^a velocities at the flux-rope footpoint inside r < and r < 2ap respectively. 
The green lines indicate the emergence velocities corresponding to the minimum shear (constant velocity, 
r < Tc) solution. Right: the blue, red and black lines correspond to the poloidal {vfjl^a ^'^'^ vertical {v1)\^^ 

1 /2 ^ ^ 

components of the optimal total velocity {v'^)2a respectively at the flux-rope footpoint inside r < 2ap. 
(both) The dotted black lines are the solar escape velocity Vq — 617kms~^, the dotted green lines are the 
photospheric sound speed Cg ~ 7.2kms~^, and the dotted blue lines are the average Alfven speed (Va)2q ■ 



agreement is nontrivial as Equations (37) and (38) indicate. Exact agreement for the magnetic field profile 
Equations (9a I and (9b) corresponds to 



dUg 
dt 



Me 
dt 



(41) 



Figure [5] shows example optimal velocity and Poynting fiux profiletj^for the 2000 September 12 CME 



assuming Bg 



B, 



4 G and Bq-^ = 500 G. The dot-dashed lines correspond to the minimum shear 



(constant velocity r < Tc) solution for the vertical velocity. The crossing point between the optimal emergence 
and constant velocity solutions occurs at the asterisks. The optimal emergence velocity profile vq, constrained 
by the energy budget, exceeds 1500 km s~^ inside the current channel r < ap and exceeds lOOkms"^ over 
most of the range < r < 3ap. The Poynting flux should be compared with the radiation omittance at 
the solar surface F ~ 6.317 x 10"'^°ergs~^ cm~^ indicated by the horizontal dotted line in the right panel of 
Figure [5] All three velocity profiles produce Poynting fiuxes which exceed the radiation omittance at the 



'^The optimal solutions for the energy flux and helicity flux are proportional with Bg in the former replaced with Ag in the 
latter and Bg/Ag = 2Bgp/apB^p = 2 Bg^/ac B(^^ ~ 2/ap, The optimal velocity and helicity flux profiles are proportional to 
those shown in Figure Isl 
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Fig. 5. — Example optimal velocity and Poynting flux profiles in the photosphere for the 2000 September 
12 CME assuming Bgc — B(^c — 4 G and B^p — 500 G. The dot-dashed lines correspond to the minimum 
shear (constant velocity r < Vc) solutions for the vertical velocity. The crossing point between the optimal 
emergence and constant velocity solutions occurs at the asterisks. The radiation emittance at the solar 
surface ~ 6.317 x 10^° erg s^^ cm~^ is indicated by the horizontal dotted line. 



solar surface and the fluxes of white light flare kernels 1 — 2 x 10^°ergs~^ cm~^ ( Neidig|["l989 ). The absence 
of strong photospheric signatures associated with CMEs suggests that flux injection cannot be responsible 
for the CME eruption. 



3.3. Discussion 



I emphasize that the 2000 September 12 event is not an extreme CME. Krall et al. (2006) describe the 



2003 October 28 CME, first observed in LASCO C3 images at 11:30 UT, that requires AUe ~ 2 x 10^3 erg and 
~ 6x10^2 Gcm-2 over Ai ~ 18 minutes with Oc = 8.1x10^ cm, S[ = 3x10^° cm and B(,c = B^^ = 3.2 G. 
The timescale of approximately 15 minutes and energy requirements 2 x 10'^^ erg are in close agreement with 



the estimates used by Manchester et al. ( 2008 1 to simulate the initiation and propagation of this eruption. 



Figure 6]shows optimal velocity estimates for the 2003 October 28 CME/ICME trajectory modeled by |KralI| 
et al.| ( 2006 ) in the same format as Figure [4] This event requires extreme photospheric velocities in the 
range v ~ 2500 — 16000 km s^^ to satisfy the energy budget and v ~ 500 — 3000 km to satisfy the helicity 
budget. 

Returning to the 2000 September 12 event, Figures |4] and [5] indicate that significant photospheric 
signatures of flux injection should be detectable in or near the current-channel radius r < 2ap. For the left 
panel in Figure [4] corresponding to ideal footpoint twisting or flux injection (emergence), all of the rms 
velocities exceed the local characteristic speeds of an MHD plasma. The optimal rms emergence velocities 

1/2 

(^c)2a sxceed the solar escape velocity vq — 617kms^^ inside of r < 2ap. Thus, gravitational forces 
cannot restrain the photospheric material in this region. Even for emergence velocities of lOOkms"^, the 
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Fig. 6.— Optimal velocity estimates for the 2003 October 28, CME/ICME trajectory modeled by [Krall 
et al. ( 2006 ) in the same format as Figure |4j 



photospheric material would require 12 minutes to return to the surface by gravity alone and this is roughly 
the timescale of the 2003 October 28, CME/ICME eruption. Consequently, a characteristic of flux injection 
should be hypersonic upflows concomitant with and preceding the eruption. Optimal photospheric velocities 
imply a sustained mass transport rate of 6 x lO^^gs"^ for just the ring-shaped regions with > 617 
maintained over roughly 40 minutes and spatial scales of tt ~ 2 x 10^° cm^ in the corona (see Figure [9] 
for the scale of the twice the current channel in the corona). These photospheric flows would eject a net 
mass of lO^'* g of which is 9 orders of magnitude larger than the typical CME mass 10^^ g as estimated 
from LASCO images. Ballistic photospheric mass outflows of this magnitude should be straight-forward to 
detect using EIT and LASCO CI and C2 coronagraphs — mass outflows of this magnitude are not observed 
by these instruments. It is worth repeating that the rms optimal emergence velocities in the photosphere are 

1/2 

largely independent of -B^ap and the cutoff scale Tc and (w^)2 a is asymptotically determined by the values 
at the base of the corona B^a^ and Oc- The left panel of Figure [4] appears to indicate that the minimum 
shear velocity (green) is more efficient than either footpoint twisting (blue) or flux injection (red) since it 
exhibits lower velocities inside r < 2 Cp. However the left panel of Figure [5] shows that this is simply because 
the profiles of the optimal twisting and emergence velocities reduce to essentially zero beyond r = Up and 
r > 1.75 flp respectively whereas the minimum shear velocity remains constant out to the cutoff scale Vc- For 
the right panel in Figure|4j corresponding to combined footpoint twisting and flux injection (emergence), flux 
emergence is only efficient when Bg Bq. For Bq ^ Bg footpoint twisting is more efficient for transferring 
energy and helicity from the photosphere to the corona. 
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Fig. 7. — EfSciencies of the driving scenarios for the 2000 September 12 CME. Left: density normahzed 
kinetic energy / = X/pp in a horizontal slice and right: ratio of the kinetic energy transport across the 
photosphere to the magnetic power requirements of the CME under the flux-injection hypothesis. The sohd 
hues correspond to rms velocities estimated from the energy budget and the dashed hues correspond to 
rms velocities estimated from the helicity budget. Blue, red, green, and black lines correspond to optimal 
twisting, optimal emergence, minimum shear (constant velocity r < r^), and optimal combined twisting and 
emergence velocity profiles respectively (blue line is not shown in right panel). 



The efficiencies of the processes may be ranked either with the magnitude of the integral 

I =2npp [ drrv"^, (42a) 



\mm ^^^^^ 



\dUe /dt\ + 2 TT pp JJ"" dr rv'^vc^^ 

where F — tt pp J^" dr r vq is the kinetic energy flux through the photosphere integrated over one foot- 
points. The former metric is simply the photospheric kinetic energy which is used to optimize the velocity 
profiles, whereas the latter metric is the ratio of the power requirements of the CME to the total energy 
transported across the photosphere under the flux-injection hypothesis. 

Figure [7] shows the efficiencies of the driving scenarios for the 2000 September 12 CME based on I, the 



kinetic energy in a horizontal slice (Equation (42a)) and e, the ratio of the magnetic power requirements 



of the CME to the total energy transported across the photosphere (Equation 42b ) under the flux-injection 
hypothesis in the left and right panels respectively. The solid lines correspond to rms velocities estimated 
from the energy budget and the dashed lines correspond to rms velocities estimated from the helicity budget. 
Blue, red, green, and black lines correspond to optimal twisting, optimal emergence, minimum shear, and 
optimal combined twisting and emergence velocities respectively — the blue line is not shown in the right 



panel because the efficiency is e = 1. Based on either criterion ( 42a[ ) or (42b) the minimum shear (constant 



velocity r < rc) is the least efficient velocity profile (shown in green) for transporting magnetic energy across 
the photosphere on the timescale of the eruption. However, the most efficient velocity is different for I and e. 
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Using the photospheric kinetic energy I, the combined twisting and emergence is the most efficient velocity 
profile (shown as the black line in the left panel) . Using the ratio of the magnetic power requirements of the 
CME to the total energy transported across the photosphere e, footpoint twisting is the most efficient velocity 
profile (e = 1) because footpoint twisting does not transport mass across the photosphere — the velocity field 
is tangent to the surface. 

The left panel Figure [T] shows that, for small vertical magnetic fields B^^ < 10 G, the kinetic energy I 
for pure emergence is slightly less than for pure footpoint twisting. For B^^ > 10 G, the kinetic energy for 
footpoint twisting less than for pure emergence. Consequently, the twisting motions dominate the kinetic 
energy of the combined twisting and emergence velocity profile since twisting motions are more efficient 
for transporting energy and helicity into the corona. The combined twisting and emergence velocity profile 
might be more efficient than the right panel of Figure[7]indicates. The rms vertical velocities for this solution 
in Figure [4] do not exceed the escape velocity and significant photospheric material could eventually return 
to the surface via strong downflows along magnetic fields. 



COMPARISON AGAINST PHOTOSPHERIC OBSERVATIONS 



The vertical photospheric plasma velocities in and near the current channel of the CME footpoints must 
be large to satisfy the flux-injection hypothesis. Such extreme velocities probably would have been detected 
in previous studies of CME eruptions, but not no such observations have been reported in the literature. 
Nonetheless, comparing expected values with observations is a necessary final step to establish the likelihood 
that a theory is compatible with nature. High spatial resolution high-cadence ('^ 1 minute) line profiles would 
be ideal for examining the footpoints of CMEs during eruptions. This suggests that limited field-of-view line- 
profile data would be the best candidate data set (see Harra et al.||2007 Imada et al.||2007 for observations 



of coronal outflows during the gradual phase of a flare). However, CME footpoints are always identified 
post-facto and are usually outside the main flux concentration of the active region that is often the focus 
of high-resolution campaigns. Thus, full-disk data are required to ensure that the dynamics of both CME 
footpoints are captured. Full-disk line-proflle data are presently scarce. One possible candidate is the Naval 
Research Laboratory Skylab He II 304 A spectroheliograms with 2" spatial resolution, but these observations 
are usually at low temporal cadence relative to flare/CME dynamics, in the wrong wavelength range for 
photospheric observations, and are difficult to disambiguate — the spatial and wavelength information is 
convolved. Finally, for direct comparison with the flux-rope model, photospheric observations concomitant 



2000 


Wood et al.||1999 


KraU et al.||2001 


Chen et al. 



that both footpoints are visible in the photosphere. 



After surveying the CMEs modeled with the flux-rope model, the 2000 September 12 CME (Chen et al 



2006 Chen & Kunkel 20101 was determined to be the best candidate because: (1) the event was front 



side and associated with a Ml.O flare, (2) the filament footpoint locations identified in (Chen et al. 2006) 
were distinct from the flare ribbons simplifying the interpretation of the observed emission lines, and (3) the 
dynamics of the filament was captured by FIT, LASCO C2 and C3 coronagraphs, and KSO Ha-observations. 
The MDI instrument was in Flarewatch mode on this day, but there is a gap in high-cadence coverage during 
the flare between 07:00 UT and 15:30 UT. However, the solar oscillations investigation (SOI) aboard SOHO 
provides continuous monitoring of Doppler velocities of low- to intermediate degree I ( Scherrer et al.|[T995 |. 
The medium-Z data are the result of Gaussian spatial filtering to mitigate spatial aliasing followed by a 
reduction in resolution from 1024 x 1024 pixels at 2"pixel~^ to 192 x 192 at roughly 10"pixel~"^. This is 



12 Sep 2000 
EIT 11:36UT 




Fig. 8. — Source region of eruption observed by EIT (Fe XII 195 A) at 11:35 UT. The erupting prominence, 
denoted by "P" appears as a dark absorption feature. The prominence footpoints are indicated "Fl" and 
"F2" and their midpoint is designated "O." The time 11:36 in the figure indicates the uncorrected start time. 



After Figure 2(a) in Chen et al. (2006) 



supplemented with level 1.8.2 2" 5-minute averaged 96-minute cadence magnetograms with the most recent 
magnetic field inter-calibrations ( Tran et al.||2005 Ulrich et al.||2009 |. 



Chen et al. (2006) identified the locations Fl and F2 of the prominence footpoints in the EIT (Fe XII 
195 A) image shown in Figure [s] Using these locations, the nominal geometry of the flux-rope model for 
the base of the corona and photosphere is diagramed in Figure |9] The top shows the EIT (Fe XII 195 A) 
image at 11:35 UT and the bottom shows the MDI (Ni I 6767.8 A) at 11:15 UT differentially de-rotated to 
11:35 UT to match the time of the EIT image using Dominic Zarro's mapping package]^ The filament and 
nominal current-channel footpoints are distinct from the fiare ribbons in the EIT image. The dark curved 
adsorbtion feature in the EIT image spanning -175" to 400" in X and -375" to -675" in Y is the filament. 

The MDI level 1.8.2 magnetograms (BLDVER18=60100) incorporate the latest sensitivities from inter- 
calibrations with the Mount Wilson Observatory (Tran et al. 2005 Ulrich et al. 2009). The line-of sight 



magnetic field was corrected for geometrical effects with a factor oi n ^ assuming that the field is purely 
vertical where 

M = ^ « Jl el/Rl, (43) 



7?2 



2Rq DqJI 



■0 



D2 



where the impact parameter qi (radians) is measured from disk center, Rq is the radius of the Sun in the 
telescope (radians or arcsecs) , and Rq / Dq is the ratio of the radius of the Sun to the distance between the 



'The SolarSoft mapping software is located at 



http: //www . Imsal . com/solarsof t/gen/idl/maps/ 
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Fig. 9. — Nominal geometry of the flux-rope model for the base of the corona and photosphere. Top: EIT 
(Fe XII 195 A) at 11:35 UT. Bottom: MDI (Ni I 6767.8 A) at 11:15 UT de-rotated to 11:35 UT to match 
the time of the EIT image. The red circles correspond to the filament footpoints in Figure [8j The small 
and large white circles correspond to the extent of twice the current-channel radius in the photosphere and 
corona respectively assuming _B^p ~ 500 G and Bi^c — 4 G. Gyan and green contours correspond to average 
magnetic fields of (-500,-400,-300) and (300,400,500) GE respectively on a scale size of 1.5 x lO^^cm^. 
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Fig. 10. — Estimate for the vertical magnetic field (flux-rope toroidal magnetic field) in the photosphere. 
Asterisks connected by dashed lines are the largest magnitude MDI values adjusted for line-of-sight projection 
and averaged over circular masks with increasing areas within the large white circles of Figure |9] Solid lines 
are Bc^^ = B^c Oc/^p estimated from the flux-rope model with i?0c = 4 G and Oc = 7.5 x 10^ cm. The top 
and bottom lines correspond to the left and right footpoints of the flux rope. See the text for discussion. 



observer and Sun center. Chen et al. (2006) argue that the filament footpoints shown in red and the nominal 



current-channel footpoints are not co-located, but are related by 

S{ = Sf + 2ac. 



(44) 



This idealized geometry is reflected in Figure |9] where the red circles correspond to the filament footpoints in 
Figure |8] and the large and small white circles correspond to the extent of twice the current-channel radius 
in the corona and photosphere respectively. The left large white circle is associated with mainly positive 
magnetic field and the right large white circle is associated with mainly negative magnetic field. 

The current-channel radius in the photosphere should be determined by enforcing consistency between 
the flux-rope model photospheric extrapolation and the MDI magnetogram though flux conservation (Equa- 

4 G with a radius a, 



cm 



tion ( 10a l). The current-channel radius at the base of the corona is — 4 G with a radius Oc = 7.5 x 10 
( Chen et al.||2006 |. If all the toroidal flux in the coronal current channel mapped to a single MDI pixel at 
disk center the resulting magnetic field would exceed 30 kG. Such extreme values of magnetic field have never 
been observed in the photosphere. Consequently, the photospheric current channel will be larger than one 
MDI pixel except perhaps very near the limb. To estimate the current-channel size in the photosphere while 
minimizing discrepancies between the photospheric observations and the flux-rope model, circular masks of 
increasing scale size were computed. These masks were convolved with the MDI values adjusted for line- 
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Fig. 11. — Temporal development of the average vertical magnetic field within r < Cp ~ 6.7 x 10® cm for 
the left (L) and right (R) filament and current-channel footpoints shown in Figure [9] The vertical dashed 
line corresponds to the rise of the M-class flare. 



of-sight projection in regions bounded by both the left and right large white circles in Figure |9] At each 
scale, the maximum magnitude magnetic fields were selected from the output of the convolution and plotted 



as asterisks connected by dashed lines in Figure 10 these values correspond to the maximum magnitude 



average magnetic field at that scale size contained within the regions bounded by the large white circles in 
Figure [o] The solid lines are B^p — B^c o-l / Op estimated from the flux-rope model with = 4 G and 
flc = 7.5 X 10^ cm. The top and bottom lines correspond to the left and right footpoints of the flux rope. 
The photospheric extrapolation of flux-rope model agrees with the MDI magnetogram where the solid and 
dashed lines cross. The photospheric area most consistent with the constraints of the flux-rope model is 
1 — 2 X 10® km^ corresponding to a current-channel radius of Op ~ 6.7 x 10® cm and an average photospheric 
toroidal held of {Bc^^) ~ 500 G. 



The geometrical relationship (|44|) proposed in Chen et al. (20061 isn't entirely consistent with the 



photospheric observations in Figure [9j For the left footpoint, the main concentration of positive flux is 
associated with the fllament footpoint not the nominal current channel. Indeed the magnetic field in the 
nominal current channel ranges from -24 to -1-24 G. For the right footpoint the magnetic field pixel ranges 
in the filament footpoint and the nominal current channel are similar ranging from —665 to —42 G and 
—600 to —62 G respectively. However, the average magnetic field magnitude in the right nominal filament 
footpoint does not exceed 300 G whereas the magnitude in the nominal current channel is somewhat higher 



as indicated by the cyan contour corresponding to -300 G contained within the small white circle. Figure 11 
shows the temporal development of the average vertical magnetic field within r < Op ~ 6.7 x 10® cm for the 
left (L) and right (R) filament and current-channel footpoints shown in Figure [o] The vertical dashed line 
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corresponds to the rise of the M-class flare. Around the time of the flare, the vertical flux is strengthening 
at the left prominence footpoint at a rate of 8 x 10"'^^ Mxs~^, but the flux in the nominal left current channel 
remains near zero over the whole five day period. The average fluxes at the right prominence and current- 
channel footpoints are flat around the eruption time. The flux-imbalance between the left and right footpoint 
indicates that the connectivity between these photospheric magnetic features is not trivial. 

While gravity certainly effects the height of the heavy prominence material potentially producing a 
difference between the position of the prominence material and the center of the current channel at the 



apex of the flux rope, there is no stated reason in Chen et al. (2006) for why the flux-rope current-channel 
and prominence footpoints aren't co-located in the photosphere]^ The footpoint separation Sf is a critical 
parameter because the height i?max of the maximum acceleration for the CME scales with Sf in the flux-rope 
model Sf/2 < iJmax ^ 3/2 ( Chen et al.||2006" Chen|[2007 ) . If S{ is not consistent with observations, then 
the flux-rope model cannot be correct. 



4.1. Doppler Data Preparation 

The MDI vector- weighted Dopplergrams were analyzed following the procedures outlined in |Snodgrass| 



(1984) and Hathaway (1992). First, the motion of the observer was removed from each Dopplergram using 
VsoHO (£»!, V') = Vr (l - e'i/2) + Vw Qi sin - Qi cosi), (45) 

where the impact parameter gi (radians) is measured from disk center and the position angle tp (radians) 
measured counterclockwise from solar north are the heliocentric radial coordinates and Vr, Vw, and Vjsi are 
the SOHO satellite velocities radial outward, westward parallel to equator, and northward along the rotation 
axis, respectively, using the appropriate keywords provided with the MDI data. 

Second, the Dopplergrams are co-registered and time-averaged. The time-averaged Dopplergram was 
fitted with orthonormal disk functions to eliminate cross-talk between coefficients ( Snodgrass 



V (Bo, g, e, $) = w (6) Rq sin $ cos + Vls (^?) + ^4if (©) + H sin 9, 
where differential rotation profile 

2 

c.(e) = ^^2„TL (sine) 



(46a) 



(46b) 



is expanded in a truncated series of even orthonormal Gegenbauer polynomials ( Morse|[l953 ), the limbshift 
function is represented by 



^LS (g) = ^ i„ (1 - cos g) , 

with (see pp. 174-175 in |Smart||1977[ ) 

g = —gi + arcsin [Dq sin (gi) / Rq] , 



(46c) 



(46d) 



^Although, the radius of the current channel at the footpoint is defined at the base of the corona the footpoints separation 
fixed by the "dense subphotospheric plasma" (see pp. 457 in |Chen||1989^ or the "massive photospheric density" ( |Chen fc] 



|Garren|1993| see pp. 2320) (see also |Krall et al.|2000| [Chen et al.|2008| . This implies that St and Si should be interpreted as 



photospheric footpoint separations (see also 



Chen & Krall 



[2004 1 
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and the meridional flow represented by 



Vmf (B) = ^ Mn Mm (sin 6) S (9) (cos $ sin 6 cos Bo - cos $ sin Bq) , 



(46e) 



are expanded in truncated series of (1 — cos g) and Fourier series of latitude respectively where the function 
classes (a;) and Mn (x) were orthonormalized by the Gram-Schmidt procedure on the interval (0, 1) by 
Snodgrass ( 1984[ ). The latitude, longitude, and solar-B angle are denoted Q, <i> and Bq respectively and 
Q = aicsm {qi / Rq) is the angle measured from the center of the Sun between the sub-Earth point and a 
point on the surface of the Sun and S (8) = +1 for 6 > and S (8) = -1 for 8 < 0. 



The function ( 46a I was fitted to the time-averaged Dopplergram with the standard deviations estimated 
from the variance of each pixel from its respective time average. The results, with a weighted = 3898.0 
and 25455 degrees of freedom, are summarized in Table [T] The first four rows of entries represent the 
best estimate of each parameter 9, their standard deviations ct, their t-scores t = 9/a, and their two- 
sided significance probability in percent that the coefficients would be larger by chance — the smaller the 
probabilities, the more significant the coefficients. The fifth and sixth rows represent th e average parameter 
estimates Ji and population standard deviations tr^ from Table 2 in Snodgrass ( 1984 1 for the data taken 
between 1967 January 1 and 1984 March 5 at the Mount Wilson Solar Observatory (WSO). The seventh 
and eighth rows represent the t-scores tgji = (^9 — Jlj / ^J^g + and their two-sided significance probability 

in percent that the best estimates 9 and the population averages JL would differ more by chance — the larger 
the probabilities, the better the two results agree. All of the best estimates are within one rms standard 



deviation of the population averages reported by Snodgrass (1984 1. 



The best estimates were used to construct a model line-of-sight velocity for the solar disk which was 
subtracted from each co-registered Dopplergram and then the Dopplergrams were de-rotated to coincide 
with the time of EIT image in Figure [9] to produce the residual velocity Awlos ■ The time average of the 



residual Doppler images is shown on the top left of Figure 12 The red circles correspond to the filament 
footpoints and the large and small white circles correspond to the extent of twice the current-channel radius 
in the corona and photosphere respectively. 

Unfortunately the MDI Dopplergrams are not absolutely calibrated. The limb-shift coefficient Lq — 
131 ± 19 ms^^ represents a constant offset for the Dopplergram series. Since the coefficient is in excellent 
agreement with measurements by ( Snodgrass|1984 ), cross-calibration against WSO with a laboratory source 
could permit absolute calibration of MDI Doppler velocities. The well known pseudo redshift in active 
regions contributes the ambiguity. The top right panel of Figure [12] shows a two-dimensional histogram 
pairing the magnitude of the magnetic field and the Doppler velocity. The magnetic data were computed by 
spatially filtering the 96-minute MDI magnetograms (also de-rotated to coincide with the time of the EIT 
image) with a two-dimensional Gaussian kernel consistent with the Medium-Z program (a = 4 in Kosovichev 



et al. 1996 1 followed by subsampling to the resolution of the Medium-^ program. Each Dopplergram was 



compared to the magnetogram that minimized differences in observation times. This panel demonstrates 
that the Doppler velocities in magnetic regions are systematically redshifted with respect to the quiet sun 
with B — 0. The histogram with \B\ < IG is Gaussian with an offset of —2.6 ± 0.2 ms~^ and a standard 
deviation of cr — 65.7 ± 0.1 nis^^. Fitting a fifth degree polynomial to the scatterplot of (Avlos) versus B 
produces the line in the right panel described by 



{vLOs){B) = J2qnB", 



(47a) 
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Fig. 12. — Top left: residual full disk Doppler velocities ft-om the Medium-Z program averaged over 24 hr. Top 
right: histogram showing the correlation between the magnitude of the magnetic field and the time-averaged 
residual Doppler velocity. Bottom: magnified view of the time-averaged residual Doppler velocities in the 
region containing the filament. The red circles correspond to the filament footpoints and the large and small 
white circles correspond to the extent of twice the current-channel radius in the corona and photosphere, 
respectively. 
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where 



qi 

92 
93 
94 
95 



-2. 2± 0.1ms" 
(5.3±0.4) X 10"^.„„ „, 
(1.75 ± 0.02) X 10^3 ms"^ 
(-1.56 ±0.08) X 10-6^ls-lG~^ 
(-4.53 ±0.09) X 10-3ms"lG~^ 



(4.9 ±0.2) X lO^^^ms^^G"^ 



(47b) 
(47c) 
(47d) 
(47e) 
(47f) 
(47g) 



However, this apparent pseudo redshift in magnetic regions is actually produced by the blueshifting of 
the quiet-sun regions caused by convective motion and the brightness velocity correlation in the convection 
cellsj^ The convective blueshift is suppressed in magnetic regions causing them to appear redshifted relative 
to the quiet sun. Consequently, the Doppler shifts may only be discussed relative to the quiet-sun motions 
which dominate the coefficient Lq and serve as the zero-point Doppler velocity of At;LOS- 

Third, the p-mode oscillations are removed from the residuals by temporal filtering. [Hathaway et all] 



(2000) employ a Gaussian weighted average of 31 Dopplergrams to reduced the p-mode signal in the 2-4 



mHz frequency band where the weights are given by 



Eli 



15 



Wi 



ith 



exp 



2a2 



exp 



2a2 



2a2 



(48a) 



(48b) 



and Ati = — to is the time difference of the central Dopplergram with a = 8 minutes and /3 = 16 minutes. 



Alternatively, nonparametric uniform B-splines h,n (j?^) — EI^i Vk Bk.2m-i (t) may be employed to 



filter temporally the data. B-splines are solutions to the optimization problem 

1 2 



N 

E 



Vi 



5v 



+ x 



oo 

/* 




dt"^ 


-oo 



(49) 



where the summation is the weighted Xmj ^vf are the uncertainties, N is number of data or knots, m is the 
spline half-order, and f] = r] (A) is a vector of B-spline parameters and a function of the global regularization 



parameter A > (Woltring 1986). The variational principle (49) may be justified by Bayesian arguments 



when the data Vi is associated with certain priors (see Craven fc Wahba]|1979 and references therein). The 
parameter A controls the balance between smoothness measured by the mth derivative and fidelity measured 
by the variance with A — > corresponding to spline interpolation with Xm — A — > cx) corresponding 
to a least squares polynomial of order P„i — m, degree d — to — 1. The boundary conditions at t — ti 
and tjv are determined by d™-hjn/dt™ — 0, where m = 1,2,..., consistent with a well-posed solution to the 
variational principle. 

B-spline basis functions have many attractive properties for filtering. For example, they do not require 
uniform sampling, they adapt to local structure, and the effect of outliers is localized because the basis 



^P. Scherrcr and B. Welsch, personal communication 2009 September; see also 



Dravins 



( 1982 I and 



Bumba & Klvaiia 



( 1995 I 
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Fig. 13. — Left: the white noise filter responses for the Gaussian weighted average filter from [Hathaway] 
et al. ( 2000 1 in green, and two m = 5th order cascaded Butterworth filter with /c = 8 mHz in blue and the 



B-spline filter in green. The filter responses were constructed with 500 realizations of a 24 hr white noise 
time series with At = 1 minute whose average spectrum is shown in black. The 2-4 mHz regime is shaded 
in light blue. Right: average power spectra of co-registered and de-rotated Doppler time-series from the 
inner 10% of the solar disk where the solid black line denotes the unfiltered data exhibiting an enhanced 
p-mode spectrum in the 2-5 mHz, the solid red line is for the smoothed results after B-spline filtering, and 
the dashed black line is for the residual data (unfiltered minus smoothed). Both left and right: the spectra 
were computed using 3 hr segments with 50% overlap. 



functions are defined on compact support. For the special case of uniformly sampled, uniformly weighted, 
periodic data, the spline smoother behaves as two cascaded mth order Butterworth filters without phase 
distortion; thus B-splines may be viewed as a generalized low pass filter adapted for nonuniform weighting^ 
uneven sampling, and non-periodic boundary conditions ( Craven fc Wahba|1979 ). Butterworth filters are an 
approximation to an ideal filter because they have maximally flat frequency response in the passband and 
are monotonically decreasing roll-off in the stop-band ( Butterworth|| 1930 ) . 



The left panel of Figure [13] shows the filter responses for the Gaussian weighted average filter from 
Hathaway et al. ( 2000 1 in green, and the two m ~ 5 order cascaded Butterworth filters with fc~8 mHz in 
blue and the B-spline filter in red with t — 1/ fc = 1250 s and 



A 



.27r 



2 At 



2m 



/At = 1.62 X lO^^s^ 



(50) 



The filter responses were constructed with 500 realizations of a 24 hour white noise time series with At = 
1 minute whose average spectrum is shown in black. The 2-4 mHz regime is shaded in light-blue. The 
B-splines exhibit a flatter power spectrum in the pass-band and a steeper roll-off in the stop-band than the 
Gaussian filter. Furthermore, the B-spline filter does not exhibit any ripples in the stop-band. Finally, the 
B-spline filter can accommodate occasional missing data or uneven sampling and simultaneously interpolate 
to a different temporal grid as part of the filtering process. Indeed, the B-spline filters were used to replace 
missing data in the central regions of four MDI Dopplergrams. The right panel of Figure [13] shows the 
average power spectra of co-registered and de-rotated Doppler time series from the inner 10% of the solar 
disk where the solid black line denotes the unfiltered data exhibiting an enhanced p-mode spectrum in the 
2-5 mHz, the solid red line is for the smoothed results after B-spline filtering, and the dashed black line is for 
the residual data (unfiltered minus smoothed). The peak power of the p-mode oscillation at / ~ 3.15 niHz 



- 28 - 



is reduced by a factor of 10" 



-10 



4.2. Photospheric Doppler Signatures During the CME 



Figure 14 shows the hne-of-sight velocity (top) and the rms line-of sight velocity (bottom) averaged 



over the current channel of radius Cp ~ 6.7 x 10 cm for the filament channel (red) and current channels 
(black) on the left (solid) and right (dashed) side of Figure [9] and the soft X-ray flux observed by GOES 8 
(blue). The M-class flare associated with the 2000 September 12 CME occurs shortly after 11:30 UT. As 
discussed in the previous section, all velocities are measured relative to the quiet sun which serves as the 
zero-point. The top plot shows that the average velocity in the left and right prominence channel and in 
the right current channel is redshifted with respect to the quiet sun. The left current channel show some 
evidence of some blueshift (upflows) with respect to the zero-point in the range of -60 to -1-10 ms~^. These 
velocities are within the standard deviation of cr = 65.7 ± 0.1 ms~^ reported for the quiet sun \B\ < 1 G 
bins in top right panel Figure |12| That velocities are quite similar to quiet sun conditions is not surprising 
because this current-channel footpoint region contains weak magnetic field ranging from -24 to 24 G. The 
weak field has already been discussed as an inconsistency between the photospheric observations and the 
simplified fiux-rope geometry which should have strong magnetic field in the current channel. Although, 
the absolute Doppler velocities are not known, they are bounded. The absolute errors must be much less 
that the passband of the filter « 7kms~^ and are likely less than 100 ms~^ because the average quiet sun 
is not moving toward the Earth at several kms~^ given the agreement for the coefficient Lq for the MDI 
Doppler measurements and ( Snodgrass|1984 ). None of the observed Doppler velocities approach magnitudes 
of lO'^ kms^^ which would be consistent with the analysis in Section [s] — these upflows would be too weak 
to transport the poloidal fiux into the corona necessary to drive the CME eruption within the context of 
the fiux-rope model. These data falsify the flux injection mechanism as a driver for this CME. However, 
an interesting feature is the sharp decrease of roughly 50ms~^ peak to minimum in the redshift of the left 
prominence channel (solid red) beginning just after the GOES X-ray flux begins to rise and lasting 1 hr 
and 35 minutes bounded by the two vertical lines. The flow velocity then appears to recover to its initial 
value over the next 2^ hr. Through this dynamic change of Doppler velocity in the filament channel, the 
magnitude remains redshifted 100 ms~^ relative to the quiet sun. To determine if dynamic change in the 
Dopper velocity of filament footpoints is a common characteristic of filament eruptions producing CMEs will 
require evaluating more events. 



COMPARISON WITH PREVIOUS WORK AND CONCLUSIONS 



The flux-injection hypothesis requires a large energy transport of 2 x lO'^'^ erg across the photosphere 
on timescales of 600-1200 s. If this hypothesis is correct, then the absence of any significant photospheric 
signature of this transport is certainly surprising given that the power requirements ~ 3 x lO'^^ergs"^ 
exceed that of the typical solar flare < lO^^ergs"^. Sudol & Harvey (2005) and Fletcher & Hudson (2008) 



have observed permanent changes in longitudinal magnetograms concomitant with solar flares, and a more 



comprehensive follow-on study by Petrie & Sudol ( 2009 ) has demonstrated that these changes are associated 



with every flare. However, these changes always lag the rise of the flare and thus would lag the peak in the 



CME acceleration which is strongly correlated with the flare rise time (Shanmugaraju et al. 2003 Zhang 



et al. 2001 2004). Furthermore, the changes in the longitudinal magnetograms are largest for events near 



the limb suggesting that they are caused by a rapid change in the angle of the magnetic field through the 
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Fig. 14. — Line-of-sight velocity (top) and the rms line-of sight velocity (bottom) averaged over the current 
channel of radius Op = 6.7 x 10^ cm for the filament channel (red) and current channels (black) on the left 
(solid) and right (dashed) and the soft X-ray flux observed by GOES 8 (blue). 

photosphere rather than a true change in photospheric magnetic field strength. 

The present study has failed to find evidence of the magnetized photospheric plasma velocities required 
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for the flux-injection hypothesis to satisfy the CME energy budget on timescales of the eruption. |Chen 



et al. (20001, Chen (2001), Krall et al. (2001), Chen & Krall (2003) and Chen & Kunkel (2010) have 



attempted to address the criticisms of the flux-injection hypothesis based on energy arguments. |Chen et al] 



(2000) cite one-dimensional simulations by Huba & Chen (1996) of a step function horizontal magnetic 



field of tens of Gauss injected into the low chromosphere. The unspecified current system generating and 
maintaining the magnetic piston is assumed to lie below the simulation boundary in the photosphere — 
this is a different situation than the diagrams in Figures [l] and [2] The simulation does not include the 
photosphere (J. Huba, personal communication 2008 December) and the plasma and neutrals only attain 
vertical velocities of ~tens of m in the low chromosphere. Chen et al. ( 2000 ) infer that these simulations 
imply photospheric upfiow velocities of merely meters per second. Can these results be used to satisfy the 
power requirements of the flux-rope model? Assuming a constant horizontal magnetic field Bh ~ 30 G 
over a large area n (S'f/2)^ ~ 2 x 10^^ cm^ surrounding each current channel and a vertical plasma velocity 
of = 10ms~^, the power supplied by this one-dimensional simulation over two footpoints would be 
dUe/dt ~ Bf^ (S'f/2)^/2 = 8 x lO^^ergs"^ which is still more than 2 orders of magnitude less than the 



power requirements of the 2000 September 12 CME (Chen et al. 2006 Chen & Kunkel 2010) which was 



associated with an Ml.O fiare and more than 3 orders of magnitude less than the power requirements of the 
2003 October 28 CME which was associated with an X17 flare ( [Krall et al.]|2006| ). 



Krall et al. ( 2000 ) consider the constraints of driving the flux-rope model by footpoint twisting. They 



conclude that footpoint twisting is inefficient and cannot reproduce the characteristics of the CME trajec- 
tories event with large photospheric poloidal velocities of lOkms"^. However, this estimate corresponds to 
velocities at the base of the corona based on their values of current-channel radius, and poloidal and toroidal 
magnetic field (see Eqs (31-33) and Tables 1 and 2 inl Krall et al.|2000 ). They coupled footpoint twisting to 
the flux-rope model by relating net current and poloidal magnetic field at the edge of the current channel to 
the number of twists at the base of the corona. However, footpoint twisting merely rearranges the current 
distribution in the current channel leaving the net current constant. Footpoint twisting should be coupled to 
the flux-rope model through the change in internal inductance as outlined in Appendix [C] not by modifying 
the net current and poloidal magnetic field at the edge of the current channel. 



Krall et al. (2001) argue that the flux-injection hypothesis implies radial plasma velocities of the order 



of Vr — Ikms""'^ at the footpoints of the flux rope implied by the inductive toroidal electric field ~ 
—c~^ dtAif,. These velocities are already substantial and should be well within the observational capabilities of 
MDI aboard SOHO and ground based observatories. However, this electric field is not relevant to estimating 
the Poynting fiux because E^xB ~ (Bg/c) dtA^r is roughly tangent to the photosphere — the vertical 
Foynting flux is zero (see Figure [2] for the geometry) . Similarly the vertical hclicity flux E X is also zero 
for the toroidal electric field. 



Chen & Krall (2003) have argued that the photospheric velocities should be highly nonuniform with 
coherence scales of less than ~ 10^ km consistent with a high f3 ^ 1 plasma outside the current channel 
although no theoretical estimate of this coherence length has been provided. High-resolution convection 
simulations exhibit rms fluctuations in the t — 1 surface of sa 30 km which is much less than the local 
photospheric pressure scale height {^casting doubt upon these assertions. Although the region far from the 
current channel is /3 ^ 1, the region containing the toroidal field and toroidal currents is /3 ~ 1. Section [3] has 
demonstrated that a significant amount of poloidal energy is injected in the region very near and interior to 
the current channel r < 2 Op. Finally, the spatial scale of the 10" pixcP^ for MDI corresponds to ~ 7500 km. 



B. Abbctt, personnel communication 2009 July. 
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The photospheric dynamics detected at this scale concomitant with the eruption at the footpoints of the 
2000 September 12 CME are not sufficient to drive the eruption. 

Emerging horizontal magnetic fields will carry mass into the chromosphere. For velocities much less than 
the escape velocity Vq = 617kms~^ this mass will flow downward as the magnetic field at the footpoints of 



the forming loops become more vertically inclined. Indeed, Ishikawa et al. (2008) have observed a strongly 
redshifted +5kms~^ Stokes V profile (downflows) at one end of an emerging horizontal magnetic feature. 



However, the turbulence argument oflChen & Krall (2003) does not seem to appreciate the very large and 



continuous average velocities~hundreds of kms~^ necessary to satisfy the power budget of the flux-rope 
model. In the limit vgp — 0, the average vertical velocities in the current channel exceed the escape velocity 
Vq and the mass could only be restrained by magnetic forces. The expected velocities also exceed the Afven 
Va and sound speed Cg in an ideal MHD plasma. Furthermore, regions of downflow embedded in the upflow 
would imply larger upflow velocities to enhance vertical transport to account for the smaller effective upflow 
area and to balance any energy transported out of the corona through the photosphere in the downflow 
regions (unless downflow regions are assumed special and contain no horizontal field). 



In a similar vein Krall et al. (2001) and Chen & Krall (2003) have argued that filling factor / — the per- 



centage of a pixel that contains magnetized atmosphere — in the photosphere is unknown and the inherent 
nonuniformity of the magnetic field in the photosphere explains the lack of observational evidence for flux 
injection. However, the filling factor can in principle be estimated by modern magnetographs by enforcing 
consistency between the Stokes / profile and the Q (linear), U (linear), and V (circular polarization) profiles 



(Keppens & Martinez Fillet 1996). Indeed, the filling factor will be a standard data product of the next 



generation vector magnetograms produced by SDO/HMl. The non-magnetic component only contributes to 
the Stokes / profile, whereas the Doppler velocities are usually determined from the V profile providing a 
relatively unambiguous estimate for the velocity of the magnetized atmosphere. Although Stokes V spectra 
may exhibit features of both upflows and downflows simultaneously ( Bellot Rubio et ar]|2001 ), the observa- 
tions in Section [4] require precise cancellation between the upflows and downflows to produce no evidence 
of significant dynamics in the flux-rope footpoint regions necessary for consistency with the flux-injection 
hypothesis. This perfect cancellation seems improbable. 

The large photospheric velocities implied by the flux-rope model power budget would shift the spectral 
lines used to estimate the magnetic field and Doppler velocity out of the pass-band of many modern telescopes. 
For example, the MDI instrument records flltergrams around the Ni I spectral line Ni I 6767.8 A with 
a 94 mA bandpass. Under normal operation, flltergrams are made at flve tuning positions separated by 
75 mA spanning 377 mA ( Scherrer et al.|l995 ). The maximum velocity that could be measured by the MDI 
instrument is u = cAA/A = 0.150/6767.8 « 7kms~^. The absence of significant changes in the magnetic 
field indicate that a large fraction of the magnetized atmosphere is not moving with velocities that exceed 
5 kms~^. 

A simple photospheric magnetic field model for fiux-rope footpoints has been developed by extrapolating 



the fiux-rope magnetic field model of Chen & Garren (1994), Chen (19961, KraU et al. (2000), and Krafi & 



Chen (2005) into the photosphere with conservation of toroidal flux and toroidal current. This is equivalent 



to magnetic fleld models implemented in simulation studies of photospheric flux injection (Chen & Huba 



2005a|b 2006 ) . This magnetic fleld model has been used to estimate the minimum photospheric velocities 
necessary to satisfy the power budget of the flux- rope model for CMEs fitted by Chen et al. (2006), Krall 



et al. (2006), and (Chen & Kunkel 2010). The flux-rope power budget requires large average poloidal or 
vertical velocities of the order of hundreds to thousands of kms~^ over large photospheric areas of 10* km^ 
to transport the necessary poloidal magnetic field into the corona on a timescale of the eruption. While 
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Chen, Krall, and Kunkel might argue that the photospheric magnetic field implemented in this study is 
oversimplified, there is agreement that enhanced photospheric activity should be detected in the region 
near the footpoints. Indeed, Krall et al. (20011 affirm "with certainty we can state that flux injection as 



discussed above should be accompanied by increased photospheric flow activity over a large spatial area, 
near the footpoints, for a period of hours during and following a CME eruption." To address this, Doppler 
and magnetic field observations at the footpoints of the 2000 September 12 CME have been analyzed. No 
significant dynamics at the flux-rope current-channel footpoints concomitant with the CME eruption have 
been detected. The flux-injection hypothesis is incompatible with these observations. 



Although, the flux-injection hypothesis has been falsifled, the flux-rope model (Chen 1989 19961 could 



remain a useful theoretical tool for modeling and interpreting CME dynamics because there are other hy- 



potheses for forming or increasing the poloidal flux of a flux rope. For example, shearing (Mikic & Linker 



1994 



Antiochos et al.|1999 Amari et al.|2000 1 , converging flows ( Forbes fc Priest|1995 ) , or nearby emerging 



flux (Chen & Shibata 2000) may convert coronal arcade held into flux-rope flelds via rapid reconnectionj^ 
Consequently, the flux-injection hypothesis is incorrect, but the flux-rope model could correctly describe the 
dynamics of an erupting CME. Finally, the flux-rope model has brought to the forefront, the paradigm that 
CMEs are current carrying coherent magnetic structures consistent with the three-part morphology observed 
by LASCO coronagraphs and predicts the scaling law that height of the CME at maximum acceleration scales 
with the footpoint separation distance S{ ( Chen et al.||2006" ). 



While the focus of this investigation has been the flux-injection hypothesis (Chen et al. 1997 2000 



[Wood et al.|[T999l [Krall et al.|[200Ti [Chen fc Kral"I|[2003l [Chen et al.||2006l [Krall et al.|[2006l [Chen fc Kunkel] 



20101, the observational component of this study places important constraints on any hypothesis that relies 
on the photosphere for the power 10^^ — lO'^'^ erg s~^ driving a CME. These hypotheses are likely incompatible 
with the present investigation, because the small velocities observed in the photospheric cannot supply the 
necessary energy on the timescale of the main acceleration phase of the CME. In contrast, the storage-release 
mechanism, where the energy is transported across the photosphere over a period of hours or days and then 
released rapidly through reconnection, is compatible with the present study. Under storage release, the 
power supplied by reconnection is limited by the reconnection rate, in part, determined by the velocity of 
the reconnection point as it unzips the overlying arcade. 
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A. THE MINIMUM PHOTOSPHERIC VELOCITIES CONSISTENT WITH THE 

POWER BUDGET 

The minimum photospheric velocities consistent with the power requirements of the flux-injection hy- 
pothesis may be found from a constrained variational principle where 

Ii - f'drrvl (Al) 

-'0 



Falsification of one hypothesis does not imply verification of another. 
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l2 = / drr Vq , 
Jo 



or 



Iz= drr (vn + v, 



is minimized subject to the integral constraint 

dUe 



Ci = 



dt 



drrBg (r) ip' (r) , 



where tp' (r) = dri^- This leads to the functional 

n=[ij' {r)f f {r)~\Be (r) ^' (r) 

and the Euler equation 



= K = constant, 



(A2) 
(A3) 

(A4) 

(A5) 
(A6) 



where A is a Lagrange multiplier. There are two limiting cases and a third general case to consider: 
Case 1: wy ^ 0, v^^ — 0, and / = B^^ . 
Solving for ip' 



V'' 



ve 



B: 



Br 



(ki + A 



(A7a) 
(A7b) 



Physical considerations require ve when r — )• which is equivalent to Bg — )• 0. This implies ki = and 



2 

\BgB^. 



Substituting Equation (A8a| into Equation (A4) determines the Lagrange multiplier 

II2OC1 



A; 



mBl^Bj^al' 



Substituting Equations (A9) and Equation (A8b| into (Al), integrating from r = 
Equations (14a) and (14b I the rms poloidal plasma velocity inside r < Op is 



2\l/2 



Case 2: v\\ 7^ 0, we = 0, and / = Bg^ 
Solving for ij)' 



dUo 



dt 



4V70 



1 



v/437 \BlB^,\al 



(A8a) 
(A8b) 

(A9) 

flp and using 
(AlO) 



CP I 



V'' 



Bl 



(k2 + A 



Bf, 



{n2 + \- 



(Alla) 
(Allb) 
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The Lagrange multiplier is determined by substituting Equation (Alia I into Equation (A4) 



X 



6 rc (85085 B^^ K2 - 170170 Ci - 144048 B^^ ' 
221 B^^al (1155 a2 - 2998 rc2) 



(A12) 



Substituting Equations ( A12) and (Allb) into Equation (A2), differentiating with respect to K2, and solving 



determines the value of K2 = corresponding to the minimum mean squared velocity over the region 
r ~ — > Tc- Integrating Equation (A2| from ?- = ^ 2ap with Equations (14a I and (14b) produces the 
rms vertical velocity inside r < 2 Cp 



dUe 
dt 



%/12516735 |B(;p| 
2Bta^ |2998B(;pr2 - 1155B(-ca?r 



(A13) 



Case 3: w = 

Solving for "0' 



0, / = 1/^/52 +52 



V'' 



Bl + B^ 



[k3 + XBe 



V = 



{K3 + XBe) 



{0,B^,-Be 



(Al4a) 
(Al4b) 



Physical considerations require vg ^ when r — )■ which is equivalent to i?e — >■ 0. This implies K3 = and 



XBg 

2 

XBg 



{Bl + Bl) 



iO,B^,~Be 



Substituting Equation (Al5a| into Equation (A4) determines the Lagrange multiplier 

36960 Ci r2 



A = 



9240 a4 - a2 ( 23984 i?L + 14421 Bl 



Cp 



Substituting Equations (AI6) and (A15b| into Equation (A3), integrating from r = 
Equations (14a) and (14b I the rms total plasma velocity inside r < Op is 



/ 2\l/2 



dUg 



dt 



2 B^p rl \/2310 ^21674^2^ + 14421 B^c B, 



Cp 



al \Bg, B^p (23984 B^ + 14421 B^c B^p) - 9240 B^ B^,\ ' 



the poloidal plasma velocity is 



/ 2\l/2 
W2a. 



dUg 



dt 



66r2 ,/30590 



BccBIp 



al \Bg, Bqp (23984 B^^ + 14421 B^^ B^p) - 9240 ^ B^ B^. 



and the vertical plasma velocity is 



dUg 



dt 



ABcp rl V12516735 



al \B^p (23984 Bj^ + 14421 B^, B^p) rl - 9240 ai B^ B^c 



(Al5a) 
(Al5b) 

(A16) 
flp and using 

(A17a) 



(Al7b) 



(A17c) 
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B. THE MINIMUM PHOTOSPHERIC VELOCITIES CONSISTENT WITH THE 

HELICITY BUDGET 

The vector potential in the upper half plane (z > 0) for an azimuthally symmetric vertical magnetic 
field is 



2 Ji (7„ r) e' ^-< /; rfrrjp (7,, r) (r) 
q;„ a (q;„) 



A^{r,z)^ — ^ ^ e (Bl) 



where an=i.2,... = 2.40483, 5.52008, ... is the nth zero of the zeroth-order Bessel function Jq (a„) — 0. This 
potential field satisfies 

V-AR,(r,C) = 0, (B2a) 
H-AR(r,0) = 0, (B2b) 
BR(r,C) = VxAR,(r,C) = V«'R,(r,C), (B2c) 
n.BR(r,0) = B^{r). (B2d) 



Ar (r,0) 



C2 = 



aB. 



Ca 



3-1-^ + 



a? 3a-^ 



r < a, 



2 



r > a. 



= — 8 7r / dr r AeR (r) -0' (r) . 



(B3) 



(B4) 



The minimum photospheric velocities consistent with the helicity budget may be found by following the 
procedure outlined in Appendix |A] using Equation (A5) with Bg Ag^. 



n^Wirf f (r)-AA,R(r) V/(r). 



Case 1: 7^ 0, = 0, and f = B^^ 

Case 2: v\\ ^ 0, vg ^ 0, and / = Bg^ 

A = 



560 C, 



437 7ra4i34^ 



dt 



2 Sec ^2998B^;7f^TT55S^ 



2310 r2 C2 



^ep^Cp'^"? (2998r2-1155a2) 



v/1155 l^cd 



Case 3: wii = and / = l/./S^ + 



2 Bg, ^12998 B<:pr2 - 1155 B^, a^J 
18480 C2 r2 



(23984 + 14421 B^J - 9240 Bj^ 



d'^g 
dt 



B^p rl V23IO ^21674 + 14421 B^^, B<^p 
2ae |S<;p (23984^2^ + 14421 B^, B^p) - 9240 B^d 



(B5) 

(B6) 
(B7) 

(B8) 
(B9) 

(BIO) 

(Bll) 
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ROTATIONAL TRANSFORM 



Krall et al. (20001 assert that footpoint twisting modifies the net current and poloidal magnetic field 



at the edge of the current channel based on a geometrical argument. Generally, the twist of a flux rope is 
estimated from the field line equations 



dr 

Br 



rdO 
Be (r) 



TZdcj) 
B4>{r)' 



(CI) 



For a circular torus, the amount of rotation of the poloidal field about the toroidal axis during the transit along 
the flux rope above the photosphere is estimated from the rotational transform (effectively the reciprocal of 
the safety factor calculation for Tokamaks) 



A0 (r) 



2Tre 



/O 

27re 



TZBg (r) 
r B^ (r) 
11 Be jr) 
rB^ (r) 



(C2a) 
(C2b) 



The amount of twist is a function of minor radius r. To add an amount of twist Si) (r) to the flux rope, the 
equation becomes 



A0 {r)+di){r) = 2 7re 



TZBg (r) 
rB^ (r) 



(5i9 (r) = 2 TT e 



TZBe jr) 
rB^ (r) 



where Be is the initial poloidal field and 

Be (r) 



Be {r) + - B^ (r) , 



27re 7^ 

Be is the final poloidal field. This equation has a similar form to Equation (3f) in 

S-Oq a, 



Krall et al. 



(2000) 



Be 



Be 



Ba 



(C3) 



(C4) 



(C5) 



27re 7^ 

where Si^a = vr is a uniform twist. However, there are important mathematical and conceptual differences 



between Equations ( C4 ) and ( C5 ) . The former Equation ( C4 ) implies that the poloidal magnetic field profile 



changes as a result of the twist whereas the latter Equation ( C5 ) relates the coefficient of Equation ( 9a 



before and after the twist. Krall et al. (20001 substitute Equation (9a) into effectively Equation (10b) to 
produce a new current 



an c 



Be 



Be 



27re n 



0c I , 



(C6) 



where 6'do = ve^t. However, the correct relationship Equation (C4) implies that Be (a) = Be {a) and 



I(t> = I(t> = Be (a) flc c/2 because (a) = 0, i.e., there is no change in net current as a consequence of the 

twisting! The effect of twisting is to modify the current distribution in the current channel and change the 

internal inductance _ 

2!^drrBj{r) 



a-'Bl 



(C7) 
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Table 1. Fit-determined CocfBcients with =3898.0 for a Reduced Set of Disk-orthonormalized 
Functions Determined from 24 hr of MDI Vector Weighted Dopplergams on 2003 September 12. Estimates 
Jb and population standard deviations dj^ from Table 2 in Snodgrass (1984) for data taken between 1967 
January 1 and 1984 March 5 at the Mount Wilson Solar Observatory. 
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